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We make new Non-LTE calculations to deduce abundances of neon from visible- 
region echelle spectra of selected Nei lines in 7 normal stars and 20 HgMn stars. We 
find that the best strong blend- free Ne line which can be used at the lower end of the 
Tgff range is A6402, although several other potentially useful Nei lines are found in the 
red region of the spectra of these stars. The mean neon abundance in the normal stars 
(log A = 8.10) is in excellent agreement with the standard abundance of neon (8.08). 
However, in HgMn stars, neon is almost universally underabundant, ranging from 
marginal deficits of 0.1-0.3 dex to underabundances of an order of magnitude or more. 
In many cases, the lines are so weak that only upper limits can be established. The most 
extreme example found is v Her with an undcrabundance of at least 1.5 dex. These 
underabundances are qualitatively expected from radiative acceleration calculations, 
which show that Ne has a very small radiative acceleration in the photosphere and is 
expected to undergo gravitational settling if mixing processes are sufficientl y weak, and 
there is no strong stellar wind. According to the theoretical predictions of Landstreeq 



et al. (199S), the low Ne abundances place an important constraint on the intensity of 
such stellar winds, which must be less than 10~^^Mq yr~^ if they are non-turbulent. 
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1 INTRODUCTION 

HgMn stars are a subclass of chemically peculiar star oc- 
cupying the spectral region corresponding approximately 
to MK types B9-B6 (10 500-16 000 K). Due to low helium 
abundances, their spectral classes are generally placed by 
observers in the A0-B8 range. Observationally, they are 
characterized by extremely low rotational velocities, weak 
or non-detectable magnetic fields and photometric variabil- 
ity, and atmospheric deficiencies of light elements (e.g.. He, 
Al, N) coupled with enhancements of the heavy elements 
(e.g., Hg, Mn, Pt, Sr, Ga). In addition, the hea vy ele- 
ment s also have non-terrestrial isotopic abundan ces ( Smith 
1997; Bohlender, Dworetsky & Jomaron 1998). The cur- 



rently favoured mechanism for explai ning these ano malies 
is the radiative diffusion hypothesis (Michaud 197C). This 



work h p.s been advanced in the form of a parameter-free 
model ( [Michaud 19861 ). 

The quiescent atmospheres of these stars makes them 
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one of the best natural laboratories for studying the compet- 
ing p rocesses of gravitational settling and radiative levita- 
tion (Vauclair & Vauclair 1982). In the absence of disrupting 



mechanisms such as convection, rotationally-induced merid- 
ional currents, high microturbulence and magnetic fields, 
certain rare elements can reach a factor of 10^ enhance- 
ment over their standard abundance. Because of the strength 
and sharpness of normally exotic spectroscopic lines, HgMn 
stars are also us eful for constraining fundamental atomic 
data (Lanz 1995). 



Although there have been many studies of individual 
HgMn stars and of the abundances of many elements across 
a sample of HgMn stars, we have been unable to find any 
papers mentioning the abundance of Ne in HgMn stars. 



with the recent exception of a paper by Adelman & Pin- 



tado (2000), in which LTE calculations established that Nei 



line strengths in k Cnc implied overabundance of 0.64 dex, 
while an undcrabundance was found in HR7245. The He-wk 
star 3 Cen A also seemed to be overabundant in Ne relative 
to the Sun (we will look again at these results in Section pf). 
However, Ne is known to exhibit strong non-LT E effects in B 
stars, according to the original investigations by Auer fc Mi 
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halas (1973), who showed that Non-LTE (NLTE) methods 
yield Ne l lines that are nearly double the strength expected 
from LTE calculations at Toh = 15 000 K. Thus the over- 
abundances reported may well be due to neglect of NLTE 
considerations. 

The lack of Ne observations in late B stars is slightly 
surprising, because Ne is an imp ortant and interesting el- 
ement. Its standard abundance (Anders & Grevesse 1989 



Grevesse, Noels & Sauval 1996), which was deduced from 



the solar wind, nebular spectrosc opy, and stellar observa- 
tions and analyses such as those of Auer fc Mihala^ , is com- 
parable with that of C, N and O. It is interesting also be- 
cause its atomic structure resembles that of He with very 
high first ionization potential (about 22 eV). Consequently, 
all its resonance lines and the ground state photoionization 
continuum are in the Lyman continuum where the stellar 
energy flux in the photosphere of late B stars is low. One 
would then expect that radiative acceleration on Ne may 



not be enough to balance gravity and Ne should sink. In- 



deed, theoretical calculations by [Landstrcct, Dolcz fc Vau 



i 



clair (1998), who considered a nonturbulent mass loss or stel- 
lar wind in the stellar envelopes, predict that there will be: 
(1) neon underabundances if the mass loss rate is less than 
1O~^''M0 yr~^; (2) neon overabundances for mass loss in 
the range 10~ — 10~^ Mq yr~ ; or (3) normal neon abun- 
dances for mass loss over lO"^'^ Mq yr~^. Also, our own cal- 
culatio ns of radiative a ccelerations in the atmospheres (puj 
daj & Dworetsky 2000) predict a pattern of general photo- 
spheric underabundances for Ne. 

In this paper, we present an abundance analysis of visi- 
ble Ne I lines based on a full NLTE treatment of the strength 
of A6402. It is the strongest unblended line of Nel in the 
visible spectrum. We also demonstrate that, to a close order 
of magnitude, several other Nel lines tend to have similar 
NLTE enhancements and can be used if spectra showing 
A6402 are not available. We find that for most HgMn stars, 
neon turns out to be underabundant, s uggesting that the 



first scenario of 
one. 



Landstreet et al. (1998 ) is the most likely 



2048 X 2048 CCDs to maximum advantage. We used both 
the unthinned phosphor-coated Orbit CCD (Dewar 13) and 
from July 1995 the thinned Ford CCD (Dewar 6), depending 
on availability as the latter was shared with the multi-object 
spectrograph on the 3-m telescope. The spectral range for 
the observations was 3800-9000 A. Typical signal-to-noise 
(S/N) per pixel in the centres of orders ranged from 75 to 
250. The Orbit CCD is cosmetically very clean, with very 
few bad pixels or columns, while the thinned Ford CCD 
contains several column defects but offers a much higher 
detector quantum efficiency in the blue. We used the Ford 
CCD whenever it was available. With the slit settings used, 
the combination of spectrographs and CCDs gave resolu- 
tions 7? ~ 46 500. Flat fields were made using the polar axis 
quartz lamp and wavelength calibrations were obtained with 
a Th-Ar comparison. 

The echelle spectra were extract ed and calibrated using 
stand ard iraf extraction packages ([Churchill 1995; Valdei 



199C), running on UCL's Starlink node. Previous measure- 



ments (Allen 1998) showed that there were no measurable 
effects of parasitic light (residual scattered light) in the line 
profiles provided that general scattered light in the adjacent 
interorder spaces was taken as the subtracted background. 
In practice the residual scattered light was less than approxi- 
mately 1 percent; we have therefore made no corrections for 
it. Allen's method is based on a direct comparison of the 
solar spectrum (as refiected from the roof of the CAT coelo- 
stat) o bserved with the HE S, with the Kitt Peak Solar Flux 
Atlas (Kurucz ct al. 1984). As the latter was obtained us- 



ing a Fourier Transform Spectrometer, it has no measurable 
parasitic light. The Kitt Peak spectrum is convolved with a 
suitable instrumental profile to match the HES data; both 
spectra must be normalized at the same points for a valid 
comparison. The ratio of summed equivalent widths of vari- 
ous features with good adjacent continuum points, in many 
different spectral orders, provides the measure of the amount 
of parasitic light. 



2 OBSERVATIONS 



Our stella r sample is based upon that of Smith fc Dworet^ 
aky (19931), who analysed lUE data on the UV resonance 



lines of iron-peak elements in 26 HgMn, 4 superficially nor- 
mal and 10 normal stars. We observed definite detections or 
det ermined upper limits for N e I in 21 of the HgMn stars in 
the smith & Dworetsky (1993) sample and 11 of the normal 



and superficially normal group. Some of the other stars in 
the two samples were lacking data in the red region or were 
cooler than 10 000 K and we did not expect to observe any 
Nel lines. Physical parameters of the stars in this study are 
given in Table nl 

All observations were obtained with the Hamilton 



Echelle Spectrograph (HES; Vogt 1987) at Lick Observa- 



tory, fed by the 0.6-m Coude Auxilliary Telescope (CAT), 
during four runs in 1994 -1997. Furth er details of the in- 
strument can be found in Misch (1997 ). Shortly before our 
observations in 1994, some of the HES optical components 
were replaced, improving the resolution and instrumental 
profile, and making it possible to use the full field of the 



3 ABUNDANCE DETERMINATION 

3.1 Stellar parameters and stellar atmospheres 

Effective temperatures and surface gravities of programme 
stars are summarized in Table nl. I n general, the parameters 
adopted follow our previous work (Smith fc Dworetsky 1993|: 



Dworetsky, Jomaron fc Smith 1998 ; Jomaron, Dworetsky & 
Allen 19991). Seven stars_are noted as double- lined spectro- 



scopic binaries in Table H; one of these (HR 1800) is better 
described as a close visual binary in which we can see evi- 
dence of the secondary spectrum as rotationally-broadened 
features. The parameters and light ratios quoted in all seven 
cases are those adopted for the primary star. Suitable light 
ra tios for other wavelength regions were found by the use 
of Kurucz (1993) model atmosphere fluxes. The light ratio 
estimated this way for A6402 is given in Table 0, where it is 
representative of the values throughout the range AA5800- 
6700. The adopted light ratio for the visual binary HR 1800 
(p = 0.243 arcsec) is 2.45 . based on A_ffp = 0.96 mag from 
The Hipparcos Catalogue (ESA 1997). Another star, 33 Gem, 
is suspected of being double- lined but there is not yet any 
information on the orbit or light ratio ( Hubrig fc Launhardt 



Table 1. Programme stars: basic data and adopted atmospheric parameters. 



HgMn stars 3 



star HD Spectral type V^ 

(kms-i) 



Toff log 3 5 

(K) (dexcms~^) (kms"'^) 



Ref. V sini 

(kms-l) 



Ref. 









Normal and 


superficially normal stars 










TrCet 


17081 


B7V 


+15 SB 


13 250 


3.80 


0.0 


[1] 


25 


[3] 


134 Tau 


38899 


B9IV 


+18 V 


10 850 


4.10 


1.6 


[1] 


30 


[3] 


rHer 


147394 


B5IV 


-14 V? 


15 000 


3.95 


0.0 


[3] 


32 


[3] 


CDra 


155763 


B6III 


-17V? 


12 900 


3.90 


2.5: 




34 


[3] 


oLyr 


172167 


AOVa 


-14V 


9 450 


4.00 


2.0 


[4] 


24 


[3] 


HR7098 


174567 


AOVs 


- 3 


10 200 


3.55 


1.0 


[3] 


11 


[3] 


21Aql 


179761 


B8 II-III 


- 5V 


13 000 


3.50 


0.2 


[1] 


17 


[3] 


HR7338 


181470 


AOIII 


-14SBO 


10 250 


3.75 


0.5 


[3] 


3 


[3] 


j/Cap 


193432 


B9.5V 


- 2V? 


10 300 


3.90 


1.6 


[1] 


27 


[3] 


HR7878 


196426 


B8IIIp 


-23 


13 050 


3.85 


1.0: 




6 


[9] 


21 Peg 


209459 


B9.5V 


- 


10 450 


3.50 


0.5 


[3] 


4 


[3] 
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87Psc 


7374 


B8III 


-16V 


13150 


4.00 


1.5 


[3] 


21.0 


[5] 


53 Tau 


27295 


B9IV 


+12SBO 


12 000 


4.25 


0.0 


[2] 


6.5 


[5] 


AiLep 


33904 


B9 IIIpHgMn 


+28 


12 800 


3.85 


0.0 


[2] 


15.5 


[5] 


HR1800 


35548 


B9 pHgSi 


- 9V? 


11050 


3.80 


0.5 


[3] 


3.0 


[5] 


33 Gem 


49606 


B7III 


+13 


14 400 


3.85 


0.5: 




22.0 


[5] 


HR2676 


53929 


B9.5III 


+ 6V? 


14 050 


3.60 


1.0: 




25.0 


[5] 


HR2844 


58661 


B9 pHgMn 


+21 V 


13 460 


3.80 


0.5: 




27.0 


[5] 


i/Cnc 


77350 


AOpSi 


-15SBO 


10 400 


3.60 


0.1 


[6] 


13 


[6] 


ftCnc 


78316 


B8IIIpMn 


+24SB10 


13 500 


3.80 


0.0 


[2] 


7 


[5] 


HR4072 


89822 


AOpSiSr:Hg: 


- 0SB2O 


10 500 


3.95 


1.0 


[7] 


3.2 


[7] 


xLup 


141556 


B9IV 


+ 5SB20 


10 750 


4.00 


0.0 


[11] 


2.0 


[7] 


tCrB 


143807 


AOpiHg: 


-19 SB 


11000 


4.00 


0.2 


[6] 


1.0 


[7] 


?jHer 


144206 


B9III 


+ 3 


12 000 


3.80 


0.6 


[1] 


9.0 


[5] 


<f>HeT 


145389 


B9p:Mn: 


-16SB10 


11650 


4.00 


0.4 


[8] 


10.1 


[5] 


HR6997 


172044 


B8 II-IIIpHg 


-26SBO 


14 500 


3.90 


1.5 


[3] 


36.0 


[5] 


112 Her 


174933 


B9 II-IIIpHg 


-20SB2O 


13100 


4.10 


0.0: 




5.5 


[12] 


HR7143 


175640 


B9III 


-26 V? 


12100 


4.00 


1.0 


[3] 


2.0 


[5] 


HR7361 


182308 


B9 IVpHgMn 


-20 V? 


13 650 


3.55 


0.0 


[3] 


8.2 


[5] 


46Aql 


186122 


B9 IIIpHgMn 


-32 


13 000 


3.65 


0.0 


[3] 


3.0 


[5] 


HR 7664 


190229 


B9 pHgMn 


-22SB1 


13 200 


3.60 


0.8 


[8] 


8.0 


[5] 


HR 7775 


193452 


AOIH 


-18" 


10 800 


3.95 


0.0 


[3] 


0.8 


[10] 



='Hoffleit & Warren (1991) cite HR7775 as SBIO although this is a confusion with /3 Cap (HR7776). 



Notes: S pectral types and radial v elocity data are from HofHeit & Warren (1991). Values of T^g and logg 
are from Smith fc Dworetsky (1993| ), or [12] in the case of 112 Her. Values of 'V and 'V?', respectiyely, indicate 
known or suspected radial velocity variables; 'SB': spect roscopic binary CSBl' and 'SB2'. re spectively, denote 



single- and double-lined systems); 'O': published orbit (see Batten, Fletcher & MacCarthy 1989). Microturbulence 
parameters f app ended by a colon (':') are approximate and were derived solely from UV Fell lines by Smith & 
Dworetsky (199: 



& Adclman (1996) 



Dworet- 



References: [11 [Adclman fc Fuhr Cl985|). [2] |Adelman (1988a|); [3] ISmith (1992f): [4] |Gigas (1986[): [5 ^ 

sky. Jomaron fc Smith fl99S)~T6rAdclrnan (198£|); 171 |Harman (1997 iHsT Adclman (1988b FTgT ::;owlcy (19"8C|); liof 
Bohlcnder. Dworetsky fc Jomaron (1998 ); [11] Wahlgrcn, Adelman fc Robinson (1994|); [12] Ryabchikova, Zakhorova 



1993); we treat it as a single star or 'average com ponent.' We lations of jeaton (199q) , who showed that his calculations 



note that Adelman, Philip fc Adelman (1996 ) also treated 



33 Ge m as a single star, noting that the question of binarity 



were in excellent agreement (within 10 per cent) with other 
recent theoretical and labor atory data such as that of Hart- 



could njt be conclusively resolved with their data. 



metz fc Schmoranzcr (1984) for the 3s-3p transitions of in- 



3.2 Atomic data for the LTK approximation 



terest in this work, and also in excellent ag reement with the 

critica lly evaluated ^/-values as given by 

(1973| ). 



Auer fc Mihalas 



As all the Ne lines in the stars observed are either weak or, 
except for the hottest stars such as r Her, not strongly satu- 
rated, the main atomic parameter of critical importance for 
LTE calculations is the oscillator strength, given as logg/ 
in Table BI We take our oscillator strengths from the calcu- 



For the radiative damping, we assumed the classical 
damping constant Tr = 2.223 x 10''/A^ s~^ [A in /xm]. This 
is a good approximation (within a factor of 2) for these lines 
as the typical lifetime of the upper levels is about 20 ns, and 
the abundances are not sensitive to the adopted values in 
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Table 2. Binary stars: adopted stellar data and light ratios 



Star A La/Lb T^ba/^oSQa Toffs/loggs L^/Lg Ref. 

(A) (K)/(cgs) (K)/(cgs) 6402A 



HR 7338 

HR 1800 

KCnc 

HR4072 

xLup 

tCrB 

112 Her 



4481 
Hp 
5480 
4520 
4520 
4520 
4520 



3.16 
2.45 
11.5 
5.45 
3.65 
2.70 
6.20 



10 250/3.8 
11050/3.8 
13 200/3.7 
10 650/3.8 
10 650/3.9 
11000/4.0 
13100/4.1 



8 500/4.0 

9 500/4.0 
8 500/4.0 

8 800/4.2 

9 200/4.2 
9 000/4.3 
8 500/4.2 



2.72 
2.34 
10.70 
5.01 
3.35 
2.46 
5.20 



[1] 
[2] 
[3] 
[4] 
[4] 
[4] 
[5] 



Note: The entry for HR 1800 is the ratio quoted for the bro adband H„ f ilter, 
which we assiime^to he the light ratio at Hfl. References: [1] Petrie (195n| ); [2] 
|eSA a997|): [3] JRvabchikova et al. I'1998|l: [4] parman (1997| and [jomaron 
Dworetsky fc Allen (199£ ); [5] [flyabchikova. et al. (199^ 



any event. Van der Waals contributions to line broadening 
are also expected to be very small; a suitable approxima- 
tion by Warner (1967) was used. For Stark broadening, we 



adopte d the rece nt experimental results of del Val, Apari- 
cio & A[ar (1999), and used an estimate of the temperature 
scaling factor proportional to T'''^ to convert their Wm at 
18 000 K to values at 12 000 K, by multiplying by an aver- 



age factor of 0.85 (Griem 1974). One line, A5852, was not 
included in their list and we adopted the simple approxi- 
mation given in CD23 data (Kurucz 1990). In general, the 
measured values which we used are about 3 times the values 
in CD23 for lines in common. We carried out worst-case sen- 



sitivity tests by varying the iel Val et al. Stark broadening 



by a factor of 2 for the strongest lines in r Her; the largest 
effect on derived abundances was less than 0.01 dex. 



log Ai'(H) — 12.00. To investigate the consistency of the re- 
sults from the selected lines, the deviations of each line from 
the mean of all the lines, \og{A/At), were calculated for each 
star where a meaningful average could be computed. These 
are summarized in Table ra where the values of log{A/Ai,) 
represent the mean deviations from the overall LTE abun- 
dance for each line. The small mean deviations imply that 
the results for each line are broadl y consistent w ith one an- 
other and the relative (jj-values of Seaton (1995). However, 
one line (A6402) is considerably stronger than all the oth- 
ers, and is well-suited for abundance determinations in the 
largest number of stars, especially those at the low-Teff end 
of the sequence and with abundances apparently below the 
standard value. In the remainder of this paper, we shall con- 
sider only this line, although future investigators may wish 
to consider some of the other lines further. 



3.3 Equivalent widths and LTE results 

Estimated abundances for several identified Nel lines were 
determined using the exact curve-of-growth technique in the 
LTE approximation. We measured the equivalent widths, 
Wx , of Ne absorption lines in the programme sp ectra by nu- 
merical integration in the dipso v3.5 package (Howarth et 
al. 1996) and compared them to the calculated values for 



each line, whi ch were generated by our s pectrum-syn thesis 
code UCLSYN dSmith fc Dworetsky 1988[ ^mith 1992| ). The 



necessary atmospheric parameters given in Table 1 
logff, and microt urbulence (5) 



were taken from 



Dworetsky (1995 ) , except for 112 Her where we used the val- 



ues given by Ryabchikova, Zakhorova & Adelman (1996) 



In the cases of the seven binaries with double spectra, we 
adopted the light ratios cited in Section 3.1 in order to cor- 



rect for dilution effects. The equivalent widths (corrected for 
binarity where necessary) and LTE abundances for se veral 
lines are^iven in Table p|. We used the 2kms~^ grid of Ku- 



rucz (1993) models, interpolating to produce a model at the 
chosen T^s and log g of each star. 



We searched a list of Nel lines from Wiese, Smith & 



Glennon (1966) in the range AA5800 — 6800 and narrowed 



the list to include only those which were fairly strong, with- 
out evident blending problems, and not situated at the ends 
of echelle orders where the spectra are noisiest. In a few cases 
the quality of the spectra justify quoting equivalent widths 
to the nearest 0.1 mA. For each star a mean LTE abundance, 
weighted by equivalent width, was calculated on the scale 



3.4 A weak blending line? 

Although we have chosen our list of Ne lines to be as blend- 
free as possible , there is a pr edicted weak blending feature in 
the CD23 list (Kurucz 199C) adjacent to the important line 
Nel A6402.246: Fen at A6402.397. In programme stars (Ta- 
ble nT) with approximately standard or lower Fe abundances, 
this line would have a typical strength of about 0.5 mA, too 
small to affect our results in any significant way. However, in 
the few stars with enhanced Fe abundances such as 112 Her 



Smith &| ( ^mith fc Dworetsky 1993[ logA(Fe) = 8.40), the strongest 



iron-rich case, the possible blending effect would have raised 
the apparent abundance of Ne by 0.3 dex. It should be noted 
that for sharp-lined stars, the displacement of the blend is 
enough to make its existence apparent. The existence of this 
line with the gf -value listed remains to be confirmed; future 
work directed at refining the neon analysis should address 
the question of its actual strength with better observations. 



4 NLTE CALCULATIONS AND 
ABUNDANCES 

4.1 Calculations and Nel model atom 

The first full NLTE calculations of neon line strengths for 
Nel were made by |Auer fc Mihalas (19"73| ) using NLTE 
model atmospheres. Unfortunately for our purposes, their 
analysis was restricted to the hotter stars with Tctt > 



HgMn stars 5 

Table 3. Ne I equivalent widths (mA) and LTE abundances for normal and HgMn programme stars on the scale logA'^(H) = 12 





A5852.49 


A6096.16 


A6266.50 


A6382.99 


A6402.25 


A6598.95 


A6717.04 


Star 


Wx 


log A 


w^ 


log A 


w^ 


log A 


Wx 


log A 


w^ 


log A 


w^ 


log A 


w^ 


log A 












Normal and 


superficially normal stars 












TrCet 


8 


8.13 


10 


8.06 


15 


8.43 


18 


8.44 


39 


8.54 


- 


- 


11 


8.37 


134 Tau 


- 


- 


- 


- 


- 


- 


- 


- 


< 15 


<8.78 


- 


- 


- 


- 


rHer 


22 


8.52 


20 


8.25 


- 


- 


26 


8.41 


59 


8.67 


- 


- 


30 


8.76 


CDra 


13: 


8.49: 


20: 


8.59: 


- 


- 


24: 


8.73: 


30 


8.32 


- 


- 


- 


- 


aLyr 


- 


- 


- 


- 


- 


- 


- 


- 


< 10 


<9.18 


- 


- 


- 




HR7098 


- 


- 


- 


- 


- 


- 


- 


- 


1.8: 


7.55: 


- 


- 


- 


- 


21Aql 


10 


8.21 


12 


8.12 


9 


8.07 


12 


8.12 


35 


8.38 


- 


- 


13 


8.43 


HR7338" 


- 


- 


- 


- 


- 


- 


- 


- 


4.0 


8.05 


- 


- 


- 


- 


uCa.p 


- 


- 


- 


- 


- 


- 


- 


- 


< 10 


<8.69 


- 


- 


- 


- 


HR7878 


8 


8.18 


14 


8.33 


8 


8.10 


16 


8.42 


30 


8.34 


4 


7.90 


10 


8.38 


21 Peg 


5 


8.65 


4 


8.37 


< 4 


<8.52 


< 2 

HgMn 


<8.03 

stars 


9 


8.36 


< 2 


<8.40 


< 2 


<8.35 


87Psc 


- 


- 


- 


- 


- 


- 


- 


- 


< 8 


< 7.41 


- 


- 


- 


- 


53 Tau 


- 


- 


- 


- 


- 


- 


- 


- 


< 5 


< 7.63 


- 


- 


- 


- 


jiLep 


- 


- 


- 


- 


- 


- 


- 


- 


< 10 


< 7.60 


- 


- 


- 


- 


HR1800'' 


- 


- 


- 


- 


- 


- 


- 


- 


< 7 


<8.02 


- 


- 


- 


- 


33 Gem 


19 


8.47 


27 


8.52 


12 


8.08 


14 


8.04 


38 


8.22 


15 


8.34 


17 


8.41 


HR2676 


- 


- 


- 


- 


- 


- 


- 


- 


25 


7.83 


- 


- 


- 


- 


HR2844 


- 


- 


- 


- 


- 


- 


- 


- 


< 16 


< 7.71 


- 


- 


- 


- 


uCnc 


- 


- 


- 


- 


- 


- 


- 


- 


< 5 


<8.05 


- 


- 


- 


- 


KCnc" 


21 


8.69 


22 


8.53 


22 


8.65 


21 


8.49 


38 


8.44 


12 


8.38 


20 


8.71 


HR4072" 


- 


- 


- 


- 


- 


- 


- 


- 


4.1 


8.04 


- 


- 


- 


- 


xLup" 


- 


- 


- 


- 


- 


- 


- 


- 


4.1 


7.95 


- 


- 


- 


- 


tCrB" 


- 


- 


- 


- 


- 


- 


- 


- 


3.7 


7.76 


- 


- 


- 


- 


uHer 


~ 


- 


- 


- 


- 


- 


- 


- 


< 1.0 


<6.63 


- 


- 


- 


- 


(/>Her 


- 


- 


- 


- 


- 


- 


- 


- 


2.0 


7.17 


- 


- 


- 


- 


HR6997 


15 


8.30 


19 


8.24 


20: 


8.38: 


20 


8.25 


37 


8.15 


- 


- 


- 


- 


112 Her" 


2.4 


7.64 


- 


- 


1.8 


7.43 


- 


- 


7.1 


7.42 


< 1-2 


< 7.43 


< 3.6 


< 7.91 


HR7143 


- 


- 


< 1 : 


< 7.3: 


< 1 : 


< 7.4: 


<2 : 


< 7.6: 


<3.0 


< 7.21 


< 1 : 


< 7.6 : 


- 


- 


HR7361 


18 


8.47 


23 


8.44 


18 


8.38 


20 


8.33 


40 


8.35 


13 


8.30 


- 


- 


46Aql 


< 2 


< 7.4 


2: 


7.2: 


3: 


7.4: 


3: 


7.4: 


9 


7.41 


2: 


7.5: 


< 2 


< 7.4 


HR7664 


8 


8.07 


13 


8.16 


11 


8.17 


12 


8.10 


29 


8.16 


8 


8.13 


10 


8.24 


HR7775 


<3 


<8.42 
-0.02 


< 2 


<8.05 
-0.05 


< 2 


<8.22 
-0.05 


6 


8.72 
-0.04 




-0.02 


< 1 


<8.16 
-0.13 


<4 


<8.79 


log(AM.) 


+0.10 


logs/ 




-0.49 




-0.31 




-0.37 




-0.24 




-1-0.33 




-0.35 




-0.35 



'^Binaries with two spectra. The Wx values are corrected for dilution effects as described in the text. Colons (':') indicate uncertain values. 



15000 K. Recently, their calculatio ns were revis it ed an d ex- 
tended to cooler temperatures by Sigut (1999 ). jigut used 
the T — T rela tions, p article densities, and electron number 



densities from Kurucq LTE line-blanketed models and solved 
the restricted NLTE problem, i.e. only the equations of ra- 
diative tranfer and statistical equlibrium for Ne. His grid of 



equivaLmt widths also has rather large steps for our pur- 
poses: 2000 K, 0.5 dex, 0.5 dex in temperature, gravity and 
Ne abundance, respectively. 

In this section we examine in detail the temperature, 
gravity and Ne-abundance region where all our HgMn stars 
are found, solve the full NLTE problem using NLTE model 
atmospheres, and find a convenient way to represent the 
NLTE efi'ects in the Nel A6402 line so that straightforward 
interpolation via LTE models can be done. For the calcu- 
lation of NLTE atmosphere models and level populations 
we used the t lusty19 5 code described in more de tail in 
Hubeny (1988| ), and in JHubeny fc Lanz (1992| , |l995| ). Here 



Hi and Nel were treated as explicit ions, which means that 
their level populations were calculated in NLTE and their 
opacity was considered. Other elements like He, C, N, and O 
were allowed to contribute to the particle and electron num- 
ber density in LTE. Synthetic spectra and equivale nt widths 



were then calculated with the SYNSPEC42 code (Hubeny, 
Jeffery 1995|). In the following, if not stated other- 



Lanz 



wise, 'in LTE' means 'in LTE considering the NLTE model 
of the atmosphere'. 



It is not possible to list here all the input parameters 
entering the NLTE calculations and above mentioned codes. 
We will mention only those which are most crucial for this 
particular problem or are difi'erent from those which could be 
generate d interactively, e.g . by the very useful interface tool 
MODION ( Varosi et al. 1995 ), part of the tlusty package for 
creating the model of the atom from the TOPbase data. We 
provide a copy of our input model for tlusty, which can be 
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Table - 

level del 

3: ionizftion energy in cm~^, col. 4: statistical weight of the level 



Nel energy levels considered. Column 1: the Paschen Table 5. 

ignation, col. 2: the nlpqr notation of Seaton (1998), col. comparin] 



Equivalent widths of N el A6402 (mA) for two models, 
this work (D+B) with kuev & Mihalas (1973). 



Paschen 


Seaton 


Energy 


g 


2p6 IS 


2p 


174192.4 


1. 


lS5 


3s332 


40148.6 


5. 


lS4 


3s331 


39731.2 


3. 


lS3 


3sll0 


39371.8 


1. 


lii-} 








2pio 


3p311 


25932.7 


3. 


2p9 


3p353 


24533.4 


7. 


2P8 


3p352 


24366.2 


5. 


2P7 


3p331 


24068.8 


3. 


2p6 


3p332 


23874.6 


5. 


2p5 


3pl31 


23418.3 


3. 


2P4 


3pl32 


23331.9 


5. 


2P3 


3p310 


23273.0 


1. 


2p2 


3plll 


23152.0 


3. 


2pi 


3pll0 


21219.7 


1. 


2ss 


4s332 


15589.3 


5. 


2S4 


4s331 


15394.4 


3. 


2S3 


4sllO 


14810.5 


1. 


2S2 


4slll 


14655.8 


3. 


3d6 


3d310 


12680.8 


1. 


3d6 


3d311 


12666.3 


3. 


3dl 


3d374 


12600.1 


9. 


3d4 


3d373 


12598.3 


7. 


3d3 


3d332 


12583.2 


5. 


3d2 


3d331 


12553.8 


3. 


3d'/ 


3d352 


12490.8 


5. 


34 


3d353 


12489.0 


7. 


3si"' 


3dl52 


11781.8 


5. 


Ss'l' 


3dl53 


11780.3 


7. 


3s'/ 


3dl32 


11770.5 


5. 


3s; 


3dl31 


11754.8 


3. 



downloaded by anyone who wishes to repeat the calculations 
(Dworetsky & Budaj 2000). 

The spectrum of Nel is that of an inert gas where 
the LS-coupling breaks down and terms and multiplets 
do not provide an appropriate description of the atom, 
so that at least the lower terms and multiplets must be 
split into individual levels and transition s. We considered 
explicitly the first 31 levels of Nel as Auer & Mihalas 
(19731), I plus continuum, with each of the levels treated 
separately (Table m). Paschen designations and exper- 
imental energies for the levels were taken from Moore 



(1949|).| Photoionization cross-sections for the terms wer e 
taken from the TOPbase data base (Cunto et al. 1993), 



as calculated by Hibbert & Scott (1994). We fit individual 



photoionization cross-sections to about 10-15 points using 
the MODION code. It was assumed that the photoionization 
cross-section is the same for the term as for the level but 
it was scaled to the particular lev el threshold. W e used the 
calculated oscillator strengths of Seaton (1998) as before. 



For collisional e xcitation rates we used th e van Regemorter 
formula as in |Auer & Mihalas (1973[). For collisional 



ionization we used equation (5.79) of Mihalas (1978) with 
gi — 0.1. With this input data we calculated the Nel A6402 
equivalent width for the same abundance (10~*), /-value 
(0.431), microturbulence ^ (4kms~^) and simil ar H-H e 
NLTE models (Tch = 15000, 20000, log 3 = 4) as |Auer 



TcB 


15000 


20000 




LTE 


NLTE 


LTE NLTE 


D+B 

A+M 


30 

28 


46 

45 


39 74 

40 79 



Table 6. Equivalent widths of Nel A6402 (mA) for two Kurucz 
models with T^a = 12000 fc 1700 K. logo = 4, g = 5kms-i, 
comparing this work (D+B) with pigut (1999). 



Tcff 




12000 




17000 




NLTE 


NLTE/LTE 


NLTE 


NLTE/LTE 


D+B 
Sigut 


19 

18 


1.31 


89 
86 


1.81 
1.80 



Mihalas to check and compare our calculations. (For the 
hotter model, two terms of Nell plus continuum were also 
considered explicitly.) The results are listed in Table 5 and 
are in ver y close agre ement. To compare our results with 
those by ^igut (199S| ), and to check the calculations for 
lower temperat ures, we solved the similar restricted NLTE 
problem using [Kurucz {199^ CD 13 LTE line blanketed 
models (computed with ^ = 2 kms^^) with our Nel atom 
model plus Hubeny's Hi atom model (9 explicit levels plus 
continuum) and the same Ne abuiidance (1.12 x 10"'*), 
/-value (0.428), and (, (5kms^*) as Sigut. The results are 
compared in Table and are also in very good agreement. 

We found that although microturbulence can affect the 
equivalent widths due to desaturation effects, for the ob- 
served stars (0 < ^ < 2.5) it has negligible effect on the 
atmosphere model, Nel level populations, and LTE/NLTE 
equivalent width ratio, R — Wa(lte)/W^a(nlte)- For this 
reason, our task can be considerably simplified as one can 
calculate a grid of models and LTE/NLTE ratios for only one 
S, — 0. Also, the NLTE effect of varying the Ne abundance 
is quite small. While decreasing the Ne abundance from 
the standard value (8.08) by 1.0 dex reduces the equivalent 
width considerably, it raises the corresponding LTE/NLTE 
ratio R by only 0.03 ± 0.01. Consequently, the main results 
can be gathered into Table M listing _R(Tcff, logg) for stan- 
dard Ne abundance and zero microturbulence. The ratio R 
ranges from around 0.6 — 0.7 at the high Tctt end of the 
HgMn domain to nearly 1.0 at the cool end of the sequence 
(where t he lines of neon disappear and can no longer be 
studied), pworetsky fc Budaj (2000| ) provide a short For- 
tran77 code to interpolate in the grid (including the small 
effects of abundance) for anyone who wishes to undertake 
their own interpola tions for R. 

Hubeny (1981) pointed out that apart from hydrogen, 
C I and Si I are the other most important ions to be included 
explicitly in the NLTE calculations in early A stars. We have 
checked that for this particular problem, including these el- 
ements along with He may again slightly affect the equiv- 
alent width, but it does not affect R significantly. The line 
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Table 7. LTE/NLTE equivalent width ratio R (Nel A6402). 



na log 5= 3.50 


3.75 


4.00 


4.25 


11000 


0.78 


0.81 


0.84 


0.87 


12000 


0.72 


0.75 


0.79 


0.82 


13000 


0.67 


0.70 


0.74 


0.77 


14000 


0.63 


0.66 


0.69 


0.73 


15000 


0.59 


0.62 


0.65 


0.69 



Table 8. Non-LTE abundances from Nel A6402 for normal and 
HgMn programme stars on the scale logA(H) = 12.00. Pro- 
gramme stars with high abundance upper limits or highly un- 
certain measurements have been omitted. 



Star 


TcB 


logs 


Wx 


W'lte/H/'nlte 


log A 




Normal and superficially 


normal stars 




TrCet 


13250 


3.80 


39 


0.69 


8.18 


rHer 


15000 


3.95 


59 


0.64 


8.15 


CDra 


12900 


3.90 


30 


0.73 


8.07 


21Aql 


13000 


3.50 


35 


0.67 


8.01 


HR7338 


10250 


3.75 


4.0 


0.86 


8.04 


HR7878 


13050 


3.85 


30 


0.71 


8.04 


21 Peg 


10450 


3.50 


9 


0.82 


8.22 






HgMn stars 






87Psc 


13150 


4.00 


< 8 


0.76 


< 7.27 


53Tau 


12000 


4.25 


< 5 


0.83 


< 7.53 


fiLep 


12800 


3.85 


< 10 


0.75 


< 7.43 


HR 1800 


11050 


3.80 


< 7 


0.83 


< 7.90 


33 Gem 


14400 


3.85 


38 


0.67 


7.87 


HR2676 


14050 


3.60 


25 


0.66 


7.54 


HR2844 


13460 


3.80 


< 16 


0.72 


< 7.50 


uCnc 


10400 


3.60 


< 5 


0.84 


< 7.94 


/tCnc 


13500 


3.80 


38 


0.69 


8.08 


HR4072 


10500 


3.95 


4.1 


0.87 


7.96 


xLup 


10750 


4.00 


4.1 


0.87 


7.87 


tCrB 


11000 


4.00 


3.7 


0.85 


7.67 


uHer 


12000 


3.80 


< 1.0 


0.80 


<6.53 


</.Her 


11650 


4.00 


2.0 


0.83 


7.08 


HR6997 


14500 


3.90 


37 


0.67 


7.82 


112 Her 


13100 


4.10 


7.1 


0.77 


7.28 


HR7143 


12100 


4.00 


< 3 


0.81 


< 7.10 


HR7361 


13650 


3.55 


40 


0.65 


7.94 


46Aql 


13000 


3.65 


9 


0.72 


7.23 


HR7664 


13200 


3.60 


29 


0.68 


7.85 



strength results were somewhat sensitive to the colhsional 
excitation. A more detailed analysis of the c urrent preci sion 
of neon NLTE calculations can be found in pigut (1999 ) . 



4.2 Neon abundances from Non-LTE calculations 

One may obtain a 'corrected' LTE equivalent width from 
RW\{ohs) and analyse the corrected width by using a stan- 
dard LTE approach including the appropriate microturbu- 
lence in a fully line-blanketed case. We used uclsyn to cal- 
culate the abundance of Ne from equivalent widths given 
in Table |3| after scaling by R from Table [TJ The microtur- 
bulence parameters f were taken from Table |l|. The results 
are shown in Table H and plotted as a function of T^s in 
Figure |l| 

Our results yield for normal and superficially-normal 
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Figure 1. Abundances of Ne in normal stars (open circles) and 
HgMn stars (filled circles). Upper limits for abundances are in- 
dicated by arrows. The s tandard abundance f or Ne on the scale 
logA(H) = 12.00 is from [Srevesse et al. (iiji ). 



stars a mean abundance log y4(Ne) — 8.10 ± 0.03 relative to 
12.00 for H. This is in excellent agr eement with the stan - 
dard abundance of Ne, 8.08, given by Grevesse et al. (1996), 



which is essentiall y identical to the value (8 .09) in the ear- 
lier compilation of Anders & Grevesse (1989). The standard 



abundance is partly based on local galactic values (stars and 
nebulae) , and on application of a well-determined correction 
to solar wind and solar energetic particle values, as well as 
spectroscopy of solar prominences. We take this agreement 
as confirmation that our ratio method for LTE/NLTE equiv- 
alent width scaling works well in the Toff range of HgMn 
stars. 

The HgMn stars are, with only one exception (kCuc), 
deficient in Ne, although the deficits in a few cases are only 
marginal (0.1 — 0.2 dex). In many cases we are able only to 
establish upper limits for Ne abundances. The most extreme 
case is v Her, for which we have particularly good spectra 
and were able to establish an upper limit 1.5 dex below 
the solar abundance. There is no case in our sample where 
Ne has an abundance greater than the standard value. The 
results of Adelman & Pintado (2000) can also be analysed 
using our method. 



5 DISCUSSION 

The error found for the average abundance for the nor- 
mal stars is based on the scatter in the results. For com- 
parison, following procedures adopte d in previous papers 
in this series on HgMn stars (e.g., Smith & Dworetsky 



199J 



Jomaron et al. 1999) we propagate uncertainties in 



adopted estimates of the errors on each parameter as fol- 
lows: ±0.25 dex in logg, ±250K in Toff, ±0.5kms"^ in the 
microturbulence (^), and ±5 percent in W\. Propagating 
these errors through the 'corrected' LTE analysis used above 
for the NelA6402 line, using a model atmosphere at T^ff = 
13 000 K, logg = 4.0, ^=1, leads to the following representa- 
tive errors in the derived Ne abundances: ±0.10 dex (logg); 
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±0.08 dex (Tcff); ±0.01 dex (C); ±0.04 dex {Wx) at standard 
abundance {W\ = 27mA), and ±0.10 dex (logg); ±0.07 dex 
(Teff); ±0.01 dex (^); ±0.03 dex {Wx) with neon underabun- 
dant by 0.5 dex iWx = 13 mA). These are very similar 
ranges; the combined expected error for one measurement is 
±0.13 dex. The s.d. for the normal stars is ±0.08 dex. This 
difference may reflect overestimates in some of the above fac- 
tors (especially A log 3) that comprise the estimated errors, 
but the two estimates are not in serious disagreement. 

We have implicitly conducted all our analyses under the 
assumption of a homogeneous depth distribution of neon in 
the photospheres of the HgMn stars, ft now seems well- 
founded to conclude that in many HgMn stars the neon 
atoms may not be distributed with a constant fraction vs. 
optical depth, due to gravitational settling, but our results 
offer no method of distinguishing clearly between a uniform 
depletion in the line forming region and an inhomogeneous 
distribution in which the total number of absorbers is about 
the same. Given the scatter in abundance from star to star, 
line strengths alone will be inadequate to prove the point 
one way or the other. 

One should consider the question of whether or not 
we may lump together the normal and superficially normal 
stars. The 'superficially normal' stars liste d in Tables h i, H 
and H were originally described as such by Cowley (1980), 
owing to the ir relatively sharp lines, but subse quent inves- 
tigati ons by smith fc Dworetsky (199C , 1993) and Smith 
(1994 H994 |199e|, |1997D have shown that HR 7338, HR7878, 



and 21 Peg have abundances which are not distinguishable 
from nor mal star s for C, N, Cr— Ni, Mg, Al, Si, Cu, Zn, Ga 
and Hg. ICowleyl thought that HR 7878 and HR 7338 were 



normal stars with no trace of peculiarity, although he sus- 
pected that 21 Peg might be related to ea rly Am stars du e 
to a weaker than expected Sc 11 line, and Sadakane (198C) 



found that both 21 Peg and HR7338 may have 'hot Am' 
characteristics, such as mild Ba and Y enhancemen ts. One 
of thes e, 21 Peg, is listed as Hg peculiar class by R enson 
et al. ( _991) without a reference cited, but Smith's (1997) 



study of Hg II lines showed that the Hg abundance is ef- 
fectively indisting uishable from that found in normal stars. 
Landstreet (1998) also describes it as a normal star. In what 



follows we assume that 21 Peg is a normal B9.5 V star and 
that HR7338 is also normal. In this work we therefore feel 
justified in includingthese three stars with the normal stars 
in Table | and Fig. 0. 



The equivalent width measures of Adelman & Pintado 



(2000Q'or 3 Gen A^ K Gnc, and HR7245 can be used with the 
results of Section W to derive approximate NLTE abundances 
for neon. For the first two of these stars, those authors give 
only the equivalent width of A5852.49 (although their ta- 
ble headings say 5842.49, which is evidently a misprint). We 
assume, based on the discussion of Section B^, that the cor- 
rection factor R may be taken to be about the same as for 
A6402, and we further assume that we can extrapolate the 
correction factor to Tcff = 17 500 for 3 Gen A, for which we 
assume R — 0.59. This well-known peculiar He-weak star 
has some characteristics similar to HgMn stars. We find a 
near-standard Ne abundance of 8.17, while k Gnc has a slight 
overabundance (8.32), and the HgMn star HR7245, which 
has a measured equivalent width of 8 mA, has a low abun- 
dance (7.32) similar to that of 112 Her. Given that our as- 
sumptions above could be subject to some uncertainty, at 



this stage we would not wish to conclude much more than 
that the abundance of neon seems consistent with the stan- 
dard value in 3 Gen A and k Gnc, but in the case of HR 7245 
we are probably on firm ground in assigning a very low abun- 
dance of neon. 

We explored briefly the question of whether the Ne 
abundances in HgMn stars depend on atmospheric parame- 
ters. It seems that largest anomalies (underabundances) are 
generally observed in the middle of the temperature range 
of HgMn stars (11500 < T^a < 13000 K; see Figure |I|). No 
apparent correlation with the surface gravity can be seen in 
Table pi It is not possible to draw any conclusions about de- 
pendences on rotational velocity as we have chosen to work 
with a selected sample of HgMn stars with fairly small v sini 
in order to ensure accurate abundance determinations. 



6 CONCLUSIONS 

We have measured the equivalent widths (or upper lim- 
its) of several Nel lines in the spectra of 11 normal late 
B stars and 21 HgMn stars in the same Teff range. These 
lines were selected after a search for lines which were well 
placed in echelle orders in the HES and which appeared not 
to have any significant blending features in the sample of 
stars studied. When analysed using LTE methods in fully 
line-blanketed atmospheres, there is a steady increase with 
Tcff in the apparent abundance above the standard value of 
log y4(Ne) = 8.08. It is apparent from previous studies that 
this is due to NLTE effects. We note that the strongest line 
of Ne in the red region, A6402, and the 6 other lines stud- 
ied, generally give concordant LTE abundances, suggesting 
that they are affected by NLTE effects by roughly the same 
amount. These lines may be of use in future investigations, 
provided further observations and NLTE calculations are 
made. 

We undertook a detailed NLTE analysis of A6402 by use 
of a full analysis including NLTE for H I and Ne I. We con- 
firm earlier studies by obtaining very similar results when 
similar inputs are used, and find that the ratio of the NLTE 
and LTE equivalent widths calculated for a given NLTE 
model atmosphere is a slowly varying function of Tch and 
surface gravity. This ratio only slightly depends on the ac- 
tual abundance of Ne, and is also very insensitive to the 
microturbulence, so it becomes possible to interpolate, in a 
table of the ratio R, the 'correction factor' by which an ob- 
served equivalent width must be scaled in order to produce 
the NLTE abundance from a much easier LTE analysis. 

The normal stars in our sample yield a mean logarith- 
mic Ne abundance of 8.10, well within our formal mean error 



of ±0.03 of the standard value, 8.08, given by Grevesse et 
al. (19961 ) . This gives us additional confidence that our mod- 
els, and the ratio method of using LTE calculations as an 
interpolation device, work satisfactorally for late B stars. 
The smallness of the scatter (s.d. ±0.08) for the individual 
stars suggests that the error budget in Section is rather 
conservative. 

It is clear from our results for the HgMn stars that the 
abundances of Ne range from standard abundance or slightly 
below, to extreme deficiencies of an order of magnitude or 
more. In several cases we have obtained only upper limits, 
and additional observations of very high quality would be 
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needed in order to attempt to detect the weak Ne lines in 
these stars. There is a tendency for the Ne abundance to 
be smallest in the middle of the HgMn effective tempera- 
ture range, but there is no dependence on surface gravity. 
There is not a single confirmed case in which Ne has an en- 
hanced abundance, which is strong evidence for the absence 



Cunto W.C., Mendoza C, Ochsenbein F., Zeippen C.J. 

1993, A&A, 275, L5 
del Val J.A., Aparicio J.A., Mar A 
Dworetsky M. M., Budaj J., 2000, 

Obser vatory. No. 81 



1999, ApJ, 513, 535 
Comm. Univ. London 



in HglV n stars of nonturbulent stellar winds (i.e. hydrogen 
mass loss rate must be < 10~ Mq yr~ ) that might com- 
pete with radiative atomic diffusion and produce accumu- 
lati ons of light elements in the photosphere, as suggested 
by Landstreet et al. (199^ ). That such winds might exist 
was studied by [Babel fc Michaud (199l| ), [Babel (199^ ) and 
Krticka fc KubAt (2000^ 



URL= http://www.ulo.ucl.ac.uk/ulo_comms 



/81/index.html 
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